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Abstract 

Using realistic low-energy model with parameters derived from the first-principles electronic 
structure calculation, we address the origin of the quasi-one-dimensional behavior in orthorhombic 
NaV204, consisting of the double chains of edge-sharing V06 octahedra. We argue that the 
geometrical aspect alone does not explain the experimentally observed anisotropy of electronic 
and magnetic properties of NaV204. Instead, we attribute the unique behavior of NaV204 to 
one particular type of the orbital ordering, which respects the orthorhombic Puma symmetry. 
This orbital ordering acts to divide all tig states into two types: the 'localized' ones, which are 
antisymmetric with respect to the mirror reflection y — >► — y, and the symmetric 'delocalized' ones. 
Thus, NaV204 can be classified as the double exchange system. The directional orientation of 
symmetric orbitals, which form the metallic band, appears to be sufficient to explain both quasi- 
one-dimensional character of interatomic magnetic interactions and the anisotropy of electrical 
resistivity. 

PACS numbers: 71.20.-b, 71.70.Gm, 75.30.Et 
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I. INTRODUCTION 



Double exchange (DE) is the key concept in the physics of strongly correlated sys- 
tems, characterizing the properties of itinerant electrons, traveling in the lattice of localized 
(atomic) spins. 1 

The situation becomes increasingly interesting when there are several degenerate orbitals, 
which become involved into the DE processes. Typically, each orbital is highly nonspherical 
and can assist the electron transfer only in some spacial directions, which can differ sub- 
stantially, depending on the shape of the occupied orbital. Thus, the itinerant electron can 
"choose" the orbital where to reside and, depending on it, choose a path for its traveling in 
the solid. Of course, the transport properties will crucially depend on such a choice. More- 
over, each guess for the orbital configuration has a feedback effect and involves some elements 
of the self-organization, typically through developing certain magnetic texture, which tends 
to stabilize the assumed orbital configuration, at least locally. Alternatively, each change of 
the magnetic state leads to the orbital-selective reconstruction of the electronic structure, 
which can be sufficient to stabilize the magnetic state.- All these phenomena are largely 
involved in the physics of pseudocubic perovskite manganites and predetermine the rich 
variety of their electronic and magnetic properties.- 

In this paper we propose that the DE physics can lead to a number of interesting effects 
when it brought in contact with the unusual crystal structure of recently synthesized mixed- 
valence (rf 15 ) NaV204 compound consisting of double chains of edge-sharing V06 octahedra 
(see Fig.[Q)i£ 

The properties of NaV204 are indeed very intriguing. It has orthorhombic space group 
Pnma, where there are four formula units in the primitive cell, and two nonequivalent 
vanadium types VI and V2, which are numbered as 1-4 and 5-8, respectively (Fig. [I]). 5 The 
vanadium atoms are surrounded by slightly distorted edge-sharing octahedra which form 
the zig-zag chains propagating in the 6-direction. 

The single-crystalline NaV204 remains metallic down to 40mK. Nevertheless, the low- 
temperature properties of NaV204 are characterized by the strong anisotropy of the electrical 
resistivity: p±/p\\ >20, where the symbols || and _L stand for the directions being parallel and 
perpendicular to the crystallographic 6-axis, respectively. In this sense, NaV204 is regarded 
as the quasi-one-dimensional (quasi- ID) compound. 
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FIG. 1: (color online). The schematic view on the crystal structure and the AFM arrangement of 
eight vanadium atoms in the primitive cell of NaV204. Eight vanadium atoms in the primitive cell 
are in the centers of octahedra and denoted by numbers (1-4 correspond to the type VI and 5-8 
- to the type V2). Small spheres in the corners of octahedra stand for the nonequivalent oxygen 
atoms. Na atoms are hidden for clarity. 

The magnetic susceptibility (x) °f NaV204 has a peak at around T N =140 K, indicating 
at an antiferromagnetic (AFM) transition. On the other hand, the fitting of x in terms of 
the Curie- Weiss law well above T N yields positive Weiss constant, indicating at the predom- 
inantly ferromagnetic (FM) character of interactions.- Furthermore, the anisotropy of elec- 
trical resistivity rapidly deteriorates above Tyy.- 1 ^ These results, together with the metallic 
character of conductivity along the 6-axis gave rise to a hypothesis that within each zig-zag 
double-chain, the intrachain interaction should be FM while the interchain interaction is 
AFM (Fig. [T]). 4 Thus, it seems that not only the electrical resistivity, but also the magnetic 
structure of NaV204 is highly anisotropic and these two anisotropies are coupled to each 
other. 

Later on, several attempts to clarify the magnetic structure of NaX^C^ have been un- 
dertaken in the experiments on the neutron scattering, 6 the muon-spin rotation,- and the 
nuclear magnetic resonance (NMR). 8 Nevertheless, the magnetic structure remains a matter 
of controversy. For example, earlier neutron-scattering measurements have been interpreted 
in terms of the incommensurate spin-density wave order with q = (0,0.191,0), in units of 
2n/b. - However, more resent muon-spin rotation studies favor the helical magnetic structure, 
where the magnetic moments lie in the ac plane and propagate along the b axis. This means 
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that some of the intrachain interactions should be AFM, which confronts the hypothesis 
based on the magnetization measurements.- Finally, there are also experimental proposals 
that NaV204 may exhibit several consequent AFM transitions below T N .-&1& 

The purpose of this work is to build some theoretical background for understanding these 
electronic and magnetic properties of NaX^C^. First, we construct a realistic low-energy 
model (Sec. [TTJ) , which aims to describe the electronic and magnetic properties of NaV204. 
Then, we analyze this model in the mean-field Hartree-Fock (HF) approximation (Sec. IIIII) . 
Although HF is a crude approximation for metallic systems, we will be able to argue, on a 
semi-quantitative level, that the main details of the electronic and magnetic structure are 
well anticipated from such a model analysis and can be attributed to the particular type of 
the orbital ordering in the metallic state. Once the orbital ordering is established, NaV2C>4 
has many things in common with the DE systems. Finally, a brief summary will be given 
in Sec. 13 



II. METHOD 



In our theoretical analysis, we follow the strategy, which was developed in the previous 
publications and which can be called as "realistic modeling of strongly correlated systems" 
(see Ref. 11 for a review). First, we calculate the electronic structure of NaV204 in the local- 
density- approximation (LDA), using the experimental parameters of the crystal structure, 
measured at T=300 K.-^ The obtained electronic structure is in a good agreement with 
the previous calculations.- Then, we construct the low-energy Hubbard-type model for the 
t 2g bands of NaV204, located near the Fermi level, and derive all parameters of such model 
from the first-principles electronic structure calculations. The details of the computational 
procedure can be found in Ref. I 111. In a number of cases, we have also constructed and solved 
the five-orbital model, which describes the behavior of both t 2g and e g bands of NaV204. 

The model includes three sets of parameters: the crystal field (CF), transfer integrals, 
and screened Coulomb interactions. 

The splitting of three t 2g levels by the CF has the following structure (in meV): (0, 46, 204) 
and (0, 10,264), at the inequivalent vanadium types VI and V2, respectively. Two lowest 
levels are nearly degenerate and not disrupted by the crystal distortion. This type of the 
CF splitting is favored by the DE interactions and minimizes the kinetic energy of electrons. 
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Nevertheless, two lower orbitals can still mix with the upper one by the transfer integrals. 
This mixing can be relatively strong. 

The behavior of some transfer integrals is explained in Fig. (2J in the local coordinate 
frame, which diagonalizes the CF Hamiltonian. The largest transfer integrals operate be- 
tween nearest-neighbor (NN) V- atoms in the chains 1-V (3-3') and 5-5' (8-8'). The transfer 
integrals between the chains are somewhat weaker, but still comparable with the ones within 
the chains. Thus, the structure of transfer integrals in is essentially three-dimensional, and 
alone does not explains the quasi- ID character of the electronic and magnetic properties of 
NaV 2 4 . 




/15 -9 47 \ 

17 89 43 
\180 17 131/ 



FIG. 2: (color online). Fragment of the crystal structure of NaV 2 04 with the matrices of transfer 
integrals between nearest neighbors. Each matrix is computed in the crystal field coordination 
frame, in the basis of t2 9 states. The matrix elements are measured in meV. The big yellow spheres 
stand for Na atoms. Other notations are the same as in Fig. [U 

The screened Coulomb interactions in the i 2 ^-band can be computed in two steps,— by ap- 
plying the constrained LDA and the random-phase approximation (RPA) for the screening. 
Roughly speaking, the first techniques takes into account the screening of atomic orbitals, 
while the second one - the self-screening by the same 3rf-electrons, which participate in the 
formation of other bands due to the hybridization effects. The fitting of screened interac- 
tions in terms of two Kanamori parameters yields the following characteristic values of the 
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intra-orbital Coulomb interaction IA— 3.15 (3.19) eV and the intraatomic exchange coupling 
J= 0.63 (0.63) eV, for the vanadium type VI (V2)M 

All parameters of the model Hamiltonian can be found in Ref. [L4|. After the construction, 
the model is solved in the HF approximation.— Since the HF approximation is known to 
have limitations for the metallic systems, we will also discuss possible impact of correlation 
interactions on the semi-quantitative level. 

Before going into details, we would like to emphasize that the analysis of electronic and 
magnetic properties of NaV204 can be greatly simplified by considering the symmetry as- 
pects. It appears that besides the CF splitting, which lifts the degeneracy of atomic t2 9 -levels 
and thus specifies subspace of the occupied orbitals, another important factor is the symme- 
try of orbitals, which are obtained from the diagonalization of the CF Hamiltonian at each 
V-site i. They can be schematically denoted as (0#i, 0^ ^#2)5 f° r ^ e vanadium type VI 
(i— 1-4), and {(f) l A) (j) l sil <% 2 )j f° r ^ e vanadium type V2 (i— 5-8). In the other words, these are 
the eigenvectors of the CF Hamiltonian, which corresponds to the eigenvalues (the CF split- 
ting), listed above. In these notations, (f) l s and (j) l A are symmetric and antisymmetric eigenvec- 
tors with respect to the mirror reflection y —> — y, which transform each V-site to itself (with 
the additional shift by b/2). For example, for the inequivalent V-sites 1 and 5 (see Fig. []]), 
these orbitals have the following form: <t>\i= 0.84|3z 2 — r 2 ) + 0.19|zx) + 0.51|x 2 — y 2 ) ) (j>\= 
0M\xy) - 0.33|2/2>, (j)\ 2 = 0.42|3z 2 -r 2 ) - 0.82|zx) - 0.38|x 2 -2/ 2 ), <j>\= 0.27\xy) - 0.96|yz), 
0|i= 0.01|3z 2 -r 2 )+0.16|zx)+0.99|x 2 - ? / 2 ), and <j>% 2 = 0.41|3z 2 -r 2 )+0.90|zx)-0.15|x 2 - ? / 2 ). 
The shape of these orbitals is explained in Fig. [31 Similar orbitals at other V-sites can be 
generated by applying the symmetry operations of the space group Pnma. The fact that 
two lowest levels, which are split off by the crystal field, belong to different representations 
(correspondingly, symmetric and antisymmetric one) has very important consequences and 
greatly reduces the number of possible solutions of the electronic model, which we have 
to consider. Indeed, due to the large orbital polarization, which is driven by the on-site 
Coulomb interaction the local populations of occupied orbitals (n), corresponding to the 
d 1 ' 5 configuration of V-sites, can be either 1 or 1/2. Since these orbitals belong to different 
representations and, therefore, do not mix with each other, we have to consider only four 
possible HF solutions, corresponding to the following populations of the symmetric and anti- 
symmetric orbitals of the vanadium types VI and V2: (n^ 1 , n A) n^ 2 , n\ 2 )= (1/2, 1, 1/2, 1), 
(1, 1/2, 1/2, 1), (1/2, 1, 1, 1/2), and (1, 1/2, 1, 1/2). Moreover, the antisymmetric orbitals at 
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(d) (e) (f) 




FIG. 3: (color online). Electron density distribution, corresponding to the crystal- field orbitals (j) 1 ^ 
(a), (b), and (j)\ 2 (c) at the V-site 1; and (p\ (d), <jP sl (e), and 0| 2 (f) at the V-site 5. 

each site will simply coincide with (j) l A (the only possible antisymmetric orbitals in the t 2g 
basis), while each occupied symmetric orbital is a linear combination of (j) l si and (j) l S2) which 
is obtained in the process of iterative solution of the HF equations. Finally, for each or- 
bital configuration, we should find the magnetic solution, corresponding to the total energy 
minimum. 

The behavior of NN transfer integrals in the chains can be also understood, by considering 
the symmetry arguments. Since for the edge-sharing geometry of the V06 octahedra, the 
main contribution to the NN integrals is caused by the ddcr interactions, which are possible 
only between orbitals of the \3z 2 — r 2 ) and \x 2 —y 2 ) type, it is reasonable to expect that the 
transfer integrals in the bonds 1-V (3-3') and 5-5' (8-8 7 ) will be large between symmetric 
orbitals, which have these components, and small - between antisymmetric ones, which do 
not have them due to the symmetry constraint. This is indeed roughly consistent with 
the behavior of the transfer integrals, depicted in Fig. El Thus, from the viewpoint of the 
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DE interactions in the chains, it is more favorable energetically to make the antisymmetric 
orbit als fully localized (n=l), and to construct the metallic band from the symmetric or- 
bitals {n—1/2). This constitutes the main idea behind results of calculations, which will be 
presented in the next section. 

III. RESULTS AND DISCUSSIONS 

We start with the analysis of the FM spin ordering, and attempt to stabilize different 
orbital structures in the HF approximation, as described above. By doing so, we were able to 
obtain only two orbital configurations, corresponding to (n^ 1 , n^ 1 , n\ 2 )= (1/2, 1, 1/2, 1) 
and (1,1/2, 1/2, 1). In the following, we will call them as Ol and (92, respectively. Two other 
solutions systematically converge to one of these two. The Ol configuration appears to be 
lower in energy by about 9 meV per one formula unit - the reason will become clear below. 
The behavior of one-electron densities of states, corresponding to the orbital configurations 
Ol and 02, is explained in Fig. E2 Both solutions are half-metallic (HM). As expected for 
the Ol configuration, the majority (t)-spin ^-states are almost fully occupied (the small 
weight in the unoccupied part of the spectrum is due to the mixing between states of the 
A- and 5-symmetry, which is caused by transfer integrals - see Fig. [2j). The metallic band 
is mainly formed by the Sl-states, with the small admixture of the 52-states. For the 02 
configuration, two V-sites display different behavior: at the site VI, the t-spin bands of the 
A and 5-symmetry are half- and totally-filled, respectively, while at the site V2, this filling 
is reversed. Moreover, the 02 alignment leads to a pseudo-gap at the Fermi level. 

The parameters of interatomic magnetic interactions are listed in Table [H They were 
calculated for two orbital configurations, by applying the perturbation theory expansion with 
respect to the infinitesimal spin rotations near the FM state.— This procedure corresponds 
to the local mapping of the total energy change near each equilibrium onto the spin model 
%s = ~^2i > j Jij e i e ji where is the direction of magnetic moment at the site i. The 
advantage of this procedure is that the parameters J^- depend on the electronic structure 
of the spin-orbital state, in which they are calculated, and can be regarded as the local 
probe of this state. One can clearly see that the change of the orbital state 01^02 has a 
dramatic effect on magnetic interactions in the VI chains and results in the sharp drop of 
the the J\y interaction. Instead, the interchain interaction J 18 tends to increase due to the 
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Energy (eV) Energy (eV) 

FIG. 4: (color online). Partial densities of states, corresponding to the ferromagnetic spin ordering 
with different orbital configurations: 01 (a) and 02 (b). VI and V2 denote two inequivalent 
vanadium sites, and A, SI and S2 - contributions of one antisymmetric and two symmetric orbitals, 
obtained from the diagonalization of the crystal field. The Fermi level is at zero energy (shown by 
dashed-dotted line). 

orbital reconstruction. On the other hand, the orbital configuration of V2 does not change so 
much. Therefore, two solutions yield similar values of interatomic magnetic interactions J 55 / 
and J 56 . As a test, we have also constructed the five-orbital model (t2 9 +e g ) and performed 
similar calculations of interatomic magnetic interactions. For 01, these results are also 
shown in Table 0J One can see that three- and five-orbital models provide similar values of 
Jij. Thus, we have confirmed that the three-orbital model captures basic magnetic properties 
of NaV204 pretty well, and in the following we will focus on the analysis of only this model. 
We have also confirmed that the change of the spin ordering for a given orbital ordering 
does not significantly modify the behavior of interatomic magnetic interactions. Therefore, 
in the following we will focus on the behavior of parameters obtained in the FM state. 

Let us concentrate on the low-energy configuration 01 and discuss the origin of the AFM 
instabilities, which were observed in the experiment. ^From the analysis of (Tabled]), one 
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TABLE I: Main parameters of interatomic magnetic interactions (in meV), calculated near the 
ferromagnetic state for two orbital configurations Ol and 02. Most of the calculations have been 
performed for the three-orbital model (the so-called model'). In addition, some test calculations 
for the Ol configuration have been also performed using the five-orbital model (the so-called L t2 g J re g 
model'). 



bond 


Ol (t 2g +e g ) 


Ol (t2g) 


02 (t 2g ) 


1-1' 


35.9 


30.1 


0.5 


1-2 


2.9 


2.0 


0.4 


1-5 


9.3 


7.7 


7.1 


1-8 


9.6 


10.8 


15.4 


1-1" 


-5.8 


-5.4 


-0.6 


5-5' 


37.8 


30.0 


26.0 


5-6 


4.1 


3.8 


4.1 


5-5" 


-3.2 


-1.0 


0.4 



can see that all NN interactions are FM. Furthermore, the FM coupling within the chains 
(Jut and J 55 /) is the strongest one, while the coupling between the chains of equivalent 
V-atoms (J12 and J 56 ), is considerably weaker (and weaker than the interactions 1-5 and 1-8 
between different double chains). 

Thus, the most probable candidate for the AFM ground state should be the one, where 
the spins between the chains of equivalent vanadium atoms are coupled antiferromagnet- 
ically and the double-chains themselves are coupled ferromagnetically (i.e., corresponding 
to the tltltllt s P m alignment at the V-atoms 1-8, depicted in Fig. []]), which seems to 
consistent with the experimental suggestion. 4 In this sense, it is right to say that the quasi- 
one-dimensional magnetic structure, realized in NaV204, is due to the Ol orbital ordering. 
Moreover, the second-neighbor interactions in the chains (Ju» and J55") are AFM. Such a 
behavior can be viewed as a precursor of the helical order, which was suggested in some 
experimental studies.- 8 Note that the AFM character of the second-neighbor interactions is 
expected for the DE systems.— Nevertheless, our calculations seem to overestimate the ten- 
dency towards the ferromagnetism. For example, the FM character of interactions Jyi and 
J 56 (Tabled]) definitely favors the FM ground state, rather than the the AFM tltltllt °ne. 
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TABLE II: Results of decomposition of nearest-neighbor magnetic interactions for the orbital 
configuration Ol in terms of double exchange (Jjj E ) and superexchange (J|- E ). All values are 
in meV. The parameter of intraatomic exchange splitting A ex was set equal to 2.8 eV in order 
to satisfy the condition = J^ E +J^ E for the largest exchange coupling Jig between different 
V-types. For other bonds, this condition is satisfied only approximately. 

bond JP E JS. E 

1-1' 45.5 -18.8 

1-2 4.2 -2.4 

1-5 17.4 -9.5 

1-8 21.6 -10.9 

5-5' 44.4 -17.4 

5-6 9.1 -5.7 

The same tendency was found in the total energy calculations, where the AFM solution was 
systematically higher than the FM one (AE= 9 and 11 meV per one formula unit in the 
three- and five-orbital model, respectively). Moreover, since FM J n / (J55/) largely exceeds 
the AFM Jq// (J55"), the helical ordering cannot be realized either. 

We attribute these discrepancies to the HF approximation, which has some limitations for 
metallic systems. Below, we qualitatively discuss the role of correlation interactions, beyond 
the HF approximation, and argue that they can indeed resolve the problem. For these 
purposes, it is convenient to decompose each NN interaction in terms of double exchange 
and superexchange (SE) contributions: J^- = J? E +J^ E . Formally, such a decomposition 
can be done for the HM electronic structure, by assuming that the splitting between the 
majority (t)- and minority (|)-spin states can be described by a single parameter A ex and 
expanding each in terms of 1/A ex . Thus, A ex has a meaning of (averaged) intraatomic 
spin splitting. Then, the zeroth-order term in this expansion corresponds to J^ E , while 
the first-order term - to J s - E . The details can be found in Ref. 



171 . The parameters of 



such an expansion are listed in Table OH We treat A ex as an adjustable parameter. Then, 
by choosing A ex « 2.8 eV, the values of interatomic magnetic interactions in Table fl] 
can be approximated reasonably well by the combination J^- « Jjj E +Jf E of DE and SE 
interactions, listed in Table ITT1 Of course, the agreement is not perfect, mainly because the 
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decomposition relies on the rigid energy splitting between the f- and |-spin states, which is 
an approximation. Nevertheless, A ex « 2.8 eV seems to be a good choice for the averaged 
splitting in the HF approximation (see Fig. Hj). As expected, all DE interactions are FM, 
while SE interactions are AFM. The large FM coupling within the chains is stabilized by the 
Ol orbital ordering, which maximizes the DE interactions (for comparison, the 02 orbital 
arrangement yields Jfy = 25.5 meV and Jfy — —26.5 meV). 

Then, it is reasonable to expect that A ex will be screened by correlation interactions. 
For example, such an effect is well known in the theory of the homogeneous electron gas.— 
Similar tendency was found in the total energy calculations for the low-energy model, which 
was constructed for the series of HM systems: the correlation interactions, treated in RPA, 
systematically decrease the effective intraatomic spin splitting and, in a number of cases, 
make this HM state unstable. 19 Thus, it is clear that the reduction of A ex , caused by the 
metallic screening effects, will strengthen J^ E , which is proportional to 1/A ex , and shift 
the balance of magnetic interactions towards the AFM coupling (of course, provided that 
the orbital configuration will not change). Particularly, in order to make the bonds 1-2 
and 5-6 antiferromagnetic (see Table [TB . the spin splitting A ex should be reduced by about 
40%. Then, if we adopt the same scaling for and note that J^y does not depend 
on A ex , the 'screened' NN interaction J\y can be estimated as J\y « Jyy+\J\y = 14.2 
meV. In this case, one can also expect the formation of the helical magnetic state in the 
chain VI, which takes place if J\y < —AJiy. Thus, the experimental magnetic ordering 
is closely related to the Ol orbital ordering: at least on the semi-quantitative level, main 
details of the experimental magnetic structure are reflected in the behavior of interatomic 
magnetic interactions. We expect that the quantitative agreement with the experimental 
data could be obtained if one goes beyond the HF approximation and consider rigorously the 
correlation interactions. This is qualitatively consistent with calculations of the correlation 
energy in RPA, by starting from one-electron eigenvalues and eigenfunctions, obtained in 
the HF approximation. This procedure yields the following values of the total (i.e., HF 
plus correlation) energy differences between AFM and FM states: AE= 6 and 2 per one 
formula unit for the three- and five-orbital model, respectively. Thus, already the 'single-shot 
calculations' of correlation energy substantially reduce AE (especially for the five-orbital 
model). We hope that a better agreement can be obtained by treating the correlation effects 
self-consistently.— 
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Finally, let us discuss the anisotropy of electrical resistivity. The dc conductivity can be 
evaluated using the Boltzmann's equation approach (in Rydberg atomic units): 20 

df(e) 



a(3 



de 



e=e z (k) 



where Q(k) is the band dispersion (/ being the band index); v t { (k) = dei(k)/dk a (p) is 
the group velocity; a(/3)= x, or z in the orthorhombic frame; r/(k) is the quasiparticle 
lifetime; and /(e) is the Fermi distribution function. For the orthorhombic Prima symmetry, 
this tensor is diagonal {a a p = (Jaadap), and the components of the resistivity tensor are given 
simply by p aa = l/cr aa . Then, for a constant r, one can easily evaluate the rations p a a/P/3(3 
between different components and compare it with the experimental data. 

Let us start our analysis with the HM FM state. Then, for the Ol orbital ordering, we 
obtain the following ratios: p xx j ' p yy = 34 and p zz j ' p yy = 19. Thus, the resistivity is highly 
anisotropic, where two components perpendicular to the chains p xx w p zz = p± are about 19- 
34 times lager than p yy = p\\ in the chain. Therefore, the experimentally observed anisotropy 
of resistivity (p±/ p\\ > 20, Ref. |4|) can be understood already from the viewpoint of the Ol 
orbital ordering, even without invoking the anisotropic tltltllt AFM arrangement. 

The orbital ordering plays a crucial role in the observed anisotropy of the resistivity 
tensor. For example, for the 02 orbital ordering, the tensor p = \\p a p\\ is nearly isotropic 
(Pxx/Pyy ~ Pzz/ Pyy ~ 1)- Similar behavior is found in the local-spin-density approximation 
for the HM FM state (p xx /p yy = 3 and p zz /p yy = 4). 

As expected,- the AFM tltltllt s Pi n alignment additionally enhances the anisotropy of 
the resistivity tensor, yielding for the Ol orbital ordering p xx /p yy = 46 and p zz /p yy = 210. 
Therefore, one may ask which ordering is more important for the anisotropy of resistivity: 
spin or orbital one? The increase of p xx j 'p yy: in comparison with the FM alignment (for the 
same orbital ordering (91), is quite modest. However, the ratio p zz /p yy changes by factor 
10. Therefore, it is tempted to conclude that the antiferromagnetism plays at least the same 
role in the anisotropy of the electrical resistivity as the orbital ordering. Nevertheless, it 
should be also remembered that the AFM order itself is caused by the orbital order (or, at 
least these two orders occur concomitantly). Thus, it is probably right to say that this is the 
orbital ordering, which has twofold effect on the anisotropy of the electrical resistivity: direct, 
which is observed already in the FM state, and indirect, via formation of the anisotropic 
AFM structure. 
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IV. SUMMARY 



On the basis of the first-principles electronic structure calculations, we have derived an 
effective low-energy model for NaV204 and employed this model for the analysis of electronic 
and magnetic properties of this material. The obtained transfer integrals are basically three- 
dimensional and by themselves do not reproduce the quasi- ID character of the electrical 
resistivity of NaX^O^ An additional important ingredient, which yields the anisotropy of 
electronic and magnetic properties should be the orbital ordering, realized in the metallic 
state of NaV204. We have argued that the symmetry arguments greatly simplify the analysis 
of the possible orbital states in NaV^O^ and the proposed orbital ordering is one of the few 
candidates, which respects the orthorhombic Prima symmetry. The AFM spin arrangement 
additionally increases the anisotropy of electrical resistivity. However, the anisotropy caused 
by the orbital ordering is comparable with the spin contribution and is formally sufficient 
for reproducing the experimental value p±/p\\ >20.- Moreover, the AFM spin ordering itself 
is driven by the orbital ordering, which increases the FM character of interactions in the 
chains and weakens the interactions between the chains. The sizable second-neighbor AFM 
interactions in the chains can be also responsible for the formation of the helical magnetic 
structure, which was proposed recently.- Nevertheless, within the HF approximation, the 
FM solution was found to be slightly lower in energy than the AFM one. In order to resolve 
this discrepancy, we provided semi-quantitative arguments and argued that the correlation 
interactions, beyond the HF approximation, should probably reverse this tendency. 
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